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Abstract 

We calculate the direct production rate of J/ ip in top quark decays. 
The color-octet J/i/j production via t — > W + bJ/ip is shown to have 
a large branching ratio of order 1.5 x 10~ 4 , which is over an order 
of magnitude higher than that of the color-singlet J/ip production 
via t — > W + bJ/ip g g or t — > W + b Xcj 9 followed by Xcj — > J/4> 7- 
This result can be used as a powerful tool to test the importance 
of the color-octet mechanism in heavy quarkonium production. 



PACS numberfs):14.40.L& 14.65.#a 



1 



Since the discovery of charmonium in 1974, there have been a lot of attempts to interpret 
the production of these new states. Among many scenarios, the color-singlet model U gains 
more success than other alternatives^ like the color-evaporation model H [[|. Based upon 
the color-singlet model, it is possible to calculate the production rates from first principles 
by standard methods^]. Indeed, the study of heavy quarkonium production may provide a 
suitable ground to precisely test quantum chromodynamics (QCD). 

However, during the past few years, it is found that the color-singlet model also has 
some defects in describing the production of heavy quarkonium. Predictions for the S- 
wave charmonium production failed to explain the new data. In the 1992-1993 run, the 
CDF detector at Fermilab Tevatron|| gave rates for prompt ip and ip' production at large 
transverse momentum which were orders of mgnitude above the lowest order perturba- 
tive calculation within the color-singlet model 0]. Even after including the fragmentation 
contributions HQ which overwhelm the former when Pt > 6Gey|10|, there is still left a 
big gap between theory and experiment. 

Recently, a new effective field theory for bound states of heavy quark and antiquark 
was provided by Bodwin, Braaten, and Lepage flTTfl in the context of non-relativistic quan- 



tum chromdynamics(NRQCD). In this approach, the interaction operators are expanded 
in powers of v , the velocity of heavy quark and antiquark in the meson rest frame, and 
a s . Contributions of different orders in v are seperated according to the "velocity scaling" 
rules. In this new framework a heavy quarkonium state H is not solely regarded as simply 
a quark-antiquark pair but rather a superposition of a series of Fock states: 

\H{nJ PC )) = O(l)|Q0( 25+1 L J ,l)) 

+ 0{v)\QQ{ 2S+ \L±l)j,M 

+ 0(v 2 )\QQC S+1 Lj,8orl)gg) 

+ (1) 

where the angular momentum of the QQ pair in each Fock state is labeled by 2S+1 Lj with a 
color configuration of either 8 or 1. The pure QQ state in color-singlet is only the leading term 
in the above expansion. Up to and including 0(v 2 ) in the Fock state expansion in describing 
J/ip(ip') production, the color-octet matrix element < Og^( 3 Si ) > [< Of ( 3 Si) >] should 
also be taken into consideration. Although these color-octet matrix elements are suppressed 
by order of v 4 relative to the corresponding color-singlet matrix elements < 0(^( 3 Sx) > 
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[< Of ( 3 Si) >], they are enhanced by a factor of l/oi^ relative to the color-singlet process 
in the short-distance perturbative calculatuion. Therefore, the supperession in the color- 
octet matrix elements can be compensated. Treating the color-octet matrix elements as free 
paramenters, the description of high-Pr J /^(ip 1 ) production at the Tevatron can indeed be 
rescued [ Yl | |1 3| , but clearly more work is needed before the new formalism is established 



as a successful theory of quarkonium production. To this end, a list of papers have been 
published [P^][IB|P7|, and still in this paper, we suggest using another important process to 
test the color-octet quarkonium production mechanism. 

The success of the Standard Model (SM) [Tj| ]l£j [pOj suggests that the top quark must exit 
||2 1|| . Recently, from the direct search at the Tevatron, the CDF and DO groups confirmed 
the existence of a heavy top quark [gg] J23f , with a mass of (176 ±8 ±10) GeV or (199^? ±22) 
GeV. The next experimental studies will focus on the determination of its properties. Among 
others, the measurement of top quark decays into heavy quark mesons which are made 
of charm or bottom quark and antiquark, will be of special interest. In particular, the 
charmonium production in top quark decays will provide very useful information in testing 
the Standard Model. 

In the Standard Model, charmonium(e.g. J/ip) may be produced via the flavor changing 
neutral current transition t — > cg|23| followed by the charm quark fragmentation c — ► J/ipc 
or the gluon fragmentation g —> J/i/j, but the rates of these processes are very small within 
the Standard Model and hence sensitive to the new physics beyond the Standard Model 



On the other hand, however, since in the Standard Model the dominant decay mode of top 
quark is t — > W + b^2j^, the dominant direct charmonium production is expected to pro- 
ceed via t — > W + bg* with the virtual gluon g* fragmentation into charmonium. There are 
three subprocesses for J/ip production via t — > W + bg*\ (i) color-octet gluon fragmenta- 
tion g* — > J/ip\ (ii) color-siglet gluon fragmentation g* — > J/ipgg; (Hi) color-singlet gluon 
fragmentation g* — > Xcjg followed by Xcj — ► J/^J, where Xcj(J = 0, 1,2) are the P-wave 
states. 

The direct charmonium production appears at orders or over in the color-singlet 
model[27], whereas the color-octet production process given by t —>■ W + bJ/ip, as shown in 
Fig.l, is at order a 2 s . Its amplitude may be written as 

A(t^wnj/i>) = l 0^T%e« w u{P') [ 7 „ r+ ± _ mfe 7a (l - 75 ) (2) 
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+ 7 Q (l-7 5 )^r-^ iMT)M t {J/il>), (3) 

p -+- /\ — m t 

where T, P', K, and P are the 4-momenta of the top quark, b quark, W + boson, and J /if) 
respectively, and M.^(J/if>) = \A(QQ[ 3 S[ 8 ^} — > J/if>)\ is the long distance nonperturbative 
amplitude of all possible way of evolving to J/if> starting from a color-octet QQ[ 3 Si^] pair at 
short distances. It may be treated as phenomenological parameter which can be determined 



by fitting the data, e.g. from the J /if) production rate at the Tevatron|l4j]. We define 

f x = — (M 2 + 2m 2 ) [(m 2 - ml) 2 + (m 2 + m\ - 2m 2 )m 2 ] , 

/a = -2(m 2 - m 2 ) 2 (M 2 + m 2 + m 2 t ) + 2m 2 y [(2M 2 + 3m 2 )(m 2 + m\) 

- 2(M 2 + ml) 2 - Am 2 b m 2 t ] , / 3 = -m\ - m 2 - 2m 2 w , 
U = 2m 4 - Am 2 b m 2 t + 2m\ - AM 2 m 2 w + 2m 2 b m 2 w + 2m 2 t m 2 w - Ami , 
h = — (M 2 + 2m 2 ) [(m 2 - m 2 b ) 2 + (m 2 + m 2 b - 2m 2 )m 2 ] , 
/ e = 2(m? - m 2 ) 2 - 2(2M 2 - m 2 - m 2 + 2m 2 )m 2 , 
/ 7 = -2m 2 - 2m 2 , / 8 = -m 2 - m 2 - 2m 2 , 

where M = 2m c is the J /if) mass. Then, the differential decay rate for t — > W + bJ/ip is given 
by 



dx\dx 2 ' 24irM 2 mfm 2 u 

+ + hxxx 2 2 + /Vx 2 ^ + /s^i^}/ (^1^2) • ( 4 ) 



Here, the variables X\ = ml — ml — ml + 2m t E w , x 2 = 2m t Ej/^ — M 2 . The physical limits 
of Xi and x 2 are 

xf = — ^ -UM 2 - x 2 ){m 2 t + m 2 b - m 2 w - x 2 ) 

2(m| - x 2 ) 

± A^[(m 2 - x 2 ), m 2 , M 2 ]A5 {(m 2 t - x 2 ),m 2 b , m 2 ]} - M 2 , 

X2 = m 2 — (m t — M) 2 , = ml — (mt, + m w ) 2 . (5) 

Here A(x, 2) = (x — y — z) 2 — 4yz. Setting a s = 0.253, m c = l.bGeV, m^ = 4.9GeV, m t = 
176GeV@, and \M 8 (J/ip)\ 2 = 0.68 x 10" 3 GeV 2 0, we get the branching ratio of 

B(t -> H/ + 6J/^) « 1.46 x 10" 4 . (6) 

The dominant color-singlet prompt J /if) production process to be t — > W + bg* with 
g* — > J/ipgg, and g* — > Xcjfi 1 followed by Xcj — ► J/ipl, as shown in Fig. 2. We can estimate 
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the partial width following the way in Ref.[g8|. The differential decay rate of t — > W + bg* is 
similar to Eq.(|]), and can be easily obtained or found in Ref.[[29|]. With the definition 



T{g*^AX)=nt/P{g*^AX), 



(7) 



the decay distribution P(g* — > Xcjg) and P(g* — > J/ipgg) for the gluon of virtuality /x can 
be found in Ref.fll and Ref.lBOI. 



r.0L 



a 



n 2 P(g* - Xcig) 

H 2 P(g* -> Xc29) 



r(l — 3r) 5 

1 — r 
6r(l + r) 

1 — r 
2r(l + 3r + 6r 2 ; 
1 — r 

l+r • 



fi P(g* -> J/ipgg) = C s r J dxj/^ J dx 1 f(x J/ ^,x 1 ;r) 

x - 

where r = M/ /x, M is the mass of the relevant charmonium states, and 



(8) 
(9) 
(10) 

(11) 



c _8a* \R' p (0)\ 2 
p 9n M 5 ' 

The function / in Eq . (|TT|) is of the formpO 

{2 + x 2 )x 2 



f(xj/j,,xi;r) 



5a 3 s \R S (0)\ 2 
s 27tt 2 M 3 ' 



(2 + xi)xi 



'12) 



(2 - x J/lP ) 2 (l -x x - r) 2 (2 - xj/^) 2 (l -x 2 - rf 



+ 
+ 



(xj/4> 



+ 



6(1 + r - xj/^f 



(1 - x 2 - r) 2 (l - xi - r) 2 (2 - X//,/,) 2 ^(1 - x 2 - r) 2 (l - x x - r) 2 
2(l-x w )(l-r) 



, in y ( 13 ) 

(1 — x 2 — r)(l — x\ — r)r r' 

where Xj = 22^/ fi with i = J/ip, gx, g 2 are the energy fractions carried by the J/ip and two 
gluons in the g* rest frame, and then x 2 = 2 — x\ — xj/^. The limits of the X\ integration in 
Eq.flTI]) are 



(14) 



We can evaluate the total decay rate of top quark to various color-singlet charmonium states, 
A, via 



mf/4 



r(t - WZ+fy*; g* ^ AX) = J dfi 2 T(t -> W~V(/x 2 )) ■ *V - ^0- 



(15) 
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In the numerical estimation, we take a s = 0.253, m c = l.bGeV, M = 2m c , \R S 



B(t - 


+ W + b Xc0 g) 


B(xco - 


- W7) 


= 2.49 x 


B(t - 


- W + b Xcl g) 




- J/H) 


= 5.35 x KT 6 


B(t - 


- W + b X c29) 


B( X c2 - 


- J/H) 


= 1.88 x KT 6 



0.999 GeV 3 , and |^(0)| 2 = 0.125 Ge\/ 5 g, and get 

(16) 
(17) 
(18) 

B(t -> W + bJ/ipgg) = 1.39 x 10" d . (19) 

The Xcj production rates depend on the infrared cutoff. Here we take the cutoff l^ min = 2M 2 , 
which is the same as that in the fragmentation analysis^]. Adding the branching ratios 
together, we obtain the total color-singlet prompt J/ip production rate to be 8.6 x 10 -6 , 
which is about a factor of 20 smaller than that via the color-octet production mechanism. 

In conclusion, we have considered the color-octet charmonium production in top quark de- 
cays, and found the branching ratio of this dominant process t — > W + bJ/ip to be 1.46 x 10~ 4 , 
which is over an order of magnitude larger than that of color-singlet production processes. 
Such a large difference makes the process of charmonium production in top decay another 
important channel to identify color-octet qurkonium signals whenever there are enough top 
quark events at the Fermilab Tevatron, LHC (Large Hadron Collider), or NLC (Next Linear 
Collider) in the future. 



Acknowlegement 

This work was supported in part by the National Natural Science Foundation of China, 
the State Education Commission of China and the State Commission of Science and Tech- 
nology of China. 



5 



References 

[1] See: for example, G. A. Schuler, Report No. CERN-TH/95-75; Report No. CERN-TH- 
7170-94, and references therein. 

[2] T. J. Weiler, Phys. Rev. Lett. 44, 304 (1980); V. Barger, W.-Y. Kueng, and R. J. N. 
Phillips, Phys. Lett. 91B, 253 (1980), and references therein. 

[3] M. B. Einhorn and S. D. Ellis, Phys. Rev. D12, 2007 (1975); H. Fritzsch, Phys. Lett. 
B67, 217 (1977); M. Gliick, J. F. Owens, and E. Reya, Phys. Rev. D17, 2324 (1978); 
J. Babcock, D. Silvers, and S. Wolfram, Phys. Rev. D18, 162 (1978). 

[4] J. H. Kwhn, J. Kaplan, and E. G. O. Safiani, Nucl. Phys. 157, 125 (1979); B. Guberina, 
J. H. Kiihn, R. D. Peccei, and R. Ruckl, Nucl. Phys. 174, 317 (1980). 

[5] CDF Collaboration, F. Abe et a/., Phys. Rev. Lett. 69 3704(1992); Phys. Rev. Lett. 
71, 2537 (1993); Contribution to the International Conference on High Energy Physics, 
Glasgow, Scotland, 1994, Report No. FERMILAB-CONF-94/136-E ( |hep-ex/9412103; ). 

[6] E. W. N. Glover, A. D. Martin, and W. J. Stirling, Z. Phys. C38, 473 (1988). 

[7] E. Braaten et al, Phys. Lett. B333, 548 (1994). 

[8] M. Cacciari and M. Greco, Phys. Rev. Lett. 73, 1586 (1994). 

[9] D. P. Roy and K. Srindhar, Phys. Lett. B339, 141 (1994). 

[10] H. A. Doncheski, S. Fleming, and M. L. Mangano, Report No. FERMILAB-CONF- 
93/348-T. 

[11] G. T. Bodwin, E. Braaten, and G. P. Lepage, Phys. Rev. D51, 1125 (1995). 
[12] E. Braaten, and S. Fleming, Phys. Rev. Lett. 74, 3327 (1995). 
[13] M. Cacciari et al, Phys. Lett. B356, 553 (1995). 

[14] P. Cho and A. K. Leibovich, Report No. CALT-68-1988 ( |hep-ph/9505329| ) and Report 
No. CALT-68-2026 ( [hep-ph/95113T5| ). 



[15] K. Cheung , W. Y. Keung, and T. C. Yuan, Phys. Rev. Lett. 76, 877 (1996). 



6 



[16] P. Cho, Report No. CALT-68-2020, 1995. 

[17] K. W. Tang and M. Vanttinen, Report No. NORDITA-96/18p; ibid, Report No. SLAC- 
PUB-95-6931. 

[18] W. Hollik, in Proceedings of the XVI International Symposium on Lepton-Photon Inter- 
actions, Cornell University, Ithaca, N.Y., Aug. 10-15 1993; M. Swartz, in Proceedings of 
the XVI International Symposium on Lepton-Photon Interactions, Cornell University, 
Ithaca, N.Y., Aug. 10-15 1993. 

[19] G. Altarelli, in Proceedings of Interational University School of Nuclear and Partical 
Physics: Substructures of Matter as Revealed with Electroweak Probes, Schladming, 
Austria, 24 Feb - 5 Mar. 1993. 

[20] G. Altarelli, CERN-TH- 7319/94, talk at 1st International Conference on Phenomenol- 
ogy of Unification: from Present to Future, Rome, Itali, 23-26 Mar 1994. 

[21] G. L. Kane, in Proceedings of the Workshop on High Energy Phenomenology, Mexico 
City, July 1-10, 1991. 

[22] F. Abe, et al. (CDF Collaboration), Phys. Rev. Lett. 74, 2626 (1995). 

[23] S. Abachi, et al. (DO Collaboration), Phys. Rev. Lett. 74, 2632 (1995). 

[24] J.L. Diaz-Cruz et al, Phys. Rev. D41 (1990) 891; B. Dutta Roy et al, Phys. Rev. Lett. 
65 (1990) 827. 

[25] Kuang-ta Chao, Cong-feng Qiao, and Feng Yuan, in preparation. 

[26] V. Barger and Phillips, Collider Physics (Addison - Wesley, Reading, MA, 1 987). 

[27] C.-F. Qiao, C. -S. Li, and K. -T. Chao, Report No. PUTP-96-03 ( [hep-ph/ 9603275) ) . 

[28] K. Hagiwara, A. D. Martin, and W. J. Stirling, Phys. Lett. B267, 527 (1991); 316, 
631(E) (1993). 

[29] G. Couture, Phys. Rev. D40, 2927 (1989); 42, 1855(E) (1990); V. Barger, A. Stange, 
and W. -Y. Keung, Phys. Rev. D42, 1835(1990). 

[30] W.-Y. Keung, Phys. Rev. D23, 2072 (1981). 



7 



[31] E. Eichten and C. Quigg, Phys. Rev. D52, 1726 (1995). 
[32] E. Braaten and T. C. Yuan, Phys. Rev. D50, 3176 (1994). 



8 



Figure Captions 



Fig.l. Color-octet Charmonium production process in top quark decays. 

Fig. 2. Diagrams for color-singlet J/ip production, (a)via g* — > J/ipgg (b) via g* — > Xcj9 — * 
J/ip'-yg. For diagram(a) Xj = 2E i / (j, with % = J/ip, g±, g2 are the energy fractions carried by 
the decay products in g* rest frame normalized to /i = m(g*). 
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